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1.  Introduction 


Ever  since  Benoit  Mandelbrot  proposed  that  natural  objects  may  be  better 
described  with  fractal  geometry  than  with  Euclidean  geometry 
(Mandelbrot,  1983),  there  have  been  numerous  attempts  to  characterize 
objects  by  their  fractal  dimension.  An  immediate  application  of  this  idea 
arises  when  we  attempt  to  classify  objects  as  subsets  of  greyscale  images, 
whether  the  image  was  produced  by  radar,  satellite  imagery,  x-ray,  black- 
and-white  photography,  or  some  similar  means.  When  such  images  are 
digitized,  they  are  reduced  to  a  range  of  greyscale  values  (often  0  through 
255)  that  are  assigned  to  a  matrix  of  pixels.  In  this  paper,  we  propose  a 
method  by  which  an  object  may  be  characterized  by  using  a  local  fractal 
dimension  operator  on  its  digitized  greyscale  image. 

The  basis  for  this  research  comes  from  the  requirement  to  identify  man¬ 
made  objects  in  images  produced  with  synthetic  aperture  radar  (SAR) 
moxmted  on  an  aerial  platform,  and  specifically  to  screen  images  of  inter¬ 
est,  which  can  then  be  confirmed  by  more  precise  techniques.  However, 
the  procedure  may  also  be  applied  to  many  other  disciplines,  particularly 
when  numerous  or  large  images  must  be  screened  and  processing  time  is 
critical.  We  note  that  the  Lincoln  Laboratory  uses  the  fractal  dimension  of 
objects  in  binarized  SAR  imagery  to  separate  natural  from  manmade  ob¬ 
jects  (Novak,  1993). 


2.  The  Fractal  Dimension  of  Natural  Objects 

The  fractal  dimension  of  a  set  A,  as  a  subset  of  R”,  is  defined  as 


D=  lim 
e-^O 


ln[N(e) 
ln(l/£)  ' 


(1) 


where  N(e)  is  the  minimum  number  of  balls  of  diameter  e  required  to 
cover  A.  While  there  are  alternative,  and  sometimes  not  equivalent 
(Hentschel,  1983),  definitions  of  fractal  dimension,  most  are  derived  from 
examining  the  limit  of  the  ratio  of  the  log  of  a  set's  "mass"  to  the  log  of  the 
scale  of  observation,  as  the  scale  goes  to  zero.  This  definition  is  consistent 
with  the  concept  of  Euclidean  dimension,  and,  in  fact,  is  a  generalization  of 
that  concept. 

This  mathematical  definition  presents  a  theoretical  problem  when  we  are 
determining  the  fractal  dimension  of  a  natural  object.  Objects  in  nature  as 
subsets  of  (or  images  as  subjects  of  R^)  will  have  differing  fractal  di¬ 
mensions  depending  on  the  scale  of  observation.  Since  there  are  no  true 
fractals  in  nature  (Mandelbrot  (1983)  points  out  that  there  are  no  straight 
lines  or  perfect  circles  in  nature  either),  we  must  modify  the  concept  of 
fractal  dimension  if  we  are  to  apply  it  to  natural  objects. 
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This  problem  is  solved  by  the  selection  of  an  appropriate  scale,  and  for  dif¬ 
ferent  classes  of  objects,  the  range  of  "e-values"  may  vary  widely 
(Turcotte,  1992).  The  appropriate  scale  for  analyzing  greyscale  imagery  is 
predetermined  by  the  set  of  pixels  over  which  the  object  is  represented.  For 
a  greyscale  image,  e  can  never  be  taken  smaller  than  the  size  of  a  single 
pixel. 


3.  Local  Fractal  Dimension 

As  the  purpose  of  our  research  is  to  classify  and  identify  objects  within  an 
image,  we  need  the  concept  of  local  fractal  dimension.  Instead  of  comput¬ 
ing  tire  fractal  dimension  of  an  image,  we  would  like  to  compute  the  fractal 
dimension  of  subsets  of  the  image,  and  use  these  local  values  to  classify  the 
subsets.  Since  our  subset  range  has  a  lower  boxmd  of  the  size  of  one  pixel, 
we  assign  a  fractal  dimension  to  each  pixel,  which  forms  the  basis  of  our 
object  classification. 

Various  techniques  exist  for  computing  a  local  fractal  dimension,  includ¬ 
ing  mass  scaling,  power  spectral  scaling,  variance  scaling  (Moghaddam, 
1991),  computing  the  Hurst  coefficient  (Gage,  1990),  and  lacimarity  (Voss, 
1986).  While  these  methods  contain  common  elements  and  all  have  a  simi¬ 
lar  result  (the  assignment  of  a  local  fractal  dimension  to  a  single  pixel 
based  on  a  window  of  surrounding  pixels),  we  are  also  interested  in  the 
speed  of  the  algorithm  for  our  application.  This  led  us  to  select  the  mass 
scaling  technique  as  the  basis  of  our  method,  which  is  most  similar  to  the 
"box-counting"  technique  of  computing  global  fractal  dimension.  Voss 
(1986)  proposed  this  type  of  local  fractal  dimension  operator,  and  Keller 
and  Chen  (1988)  published  an  algorithm  to  implement  the  procedure. 
While  other  studies  found  this  to  be  the  least  accurate  method 
(Moghaddam,  1990),  it  was  the  only  one  that  promised  real-time  results,  a 
feature  critical  to  our  application. 


4.  Measure  of  "Roughness" 

Since  a  digitized  image  is  essentially  a  discretization  of  a  surface  in  three 
dimensions,  we  desire  an  operator  that  assigns  a  fractal  dimension  consis¬ 
tent  with  the  "roughness"  of  the  surface,  i.e.,  higher  fractal  dimensions  cor¬ 
responding  to  rougher  surfaces,  with  a  value  of  two  for  a  flat  surface  and  a 
value  of  three  for  a  space-filling  surface. 

The  goal  of  having  a  higher  fractal  dimension  assigned  to  a  "rougher"  sur¬ 
face  leads  to  an  ambiguity  with  local  fractal  dimension  operators:  While 
there  are  ways  to  measure  surface  roughness  (Milman,  1994),  what  is  lack¬ 
ing  is  a  clear  definition  of  what  constitutes  a  "rough"  surface  in  terms  of 
mass  scaling.  We  demonstrate  this  problem  by  example  with  three  5x5 
"windows"  of  pixels  about  a  central  pixel,  for  which  the  local  fractal  di¬ 
mension  is  to  be  computed. 


6 


0 

200 

0 

200 

0 

90 

110 

90 

110 

90 

50 

0 

100 

150 

200 

200 

0 

200 

0 

200 

110 

90 

110 

90 

110 

50 

0 

100 

150 

200 

0 

200 

100 

200 

0 

90 

110 

100 

110 

90 

50 

0 

100 

150 

200 

200 

0 

200 

0 

200 

110 

90 

110 

90 

110 

50 

0 

100 

150 

200 

0 

200 

0 

200 

0 

90 

110 

90 

110 

90 

50 

0 

100 

150 

200 

Case  I  Case  II  Case  III 


Figure  1.  Greyscale  values  for  three  5x5  image  segments. 

Consider  the  three  cases  shown  in  figure  1.  Case  I  is  clearly  more  rough 
than  case  II,  for  although  these  two  image  segments  have  the  same 
"shape"  about  the  central  pixel,  the  dynamic  range  is  much  greater  in  case 
I.  Therefore,  we  would  expect  any  local  fractal-dimension  operator  to  as¬ 
sign  a  higher  value  to  the  central  pixel  in  case  I. 

Similarly,  case  I  is  rougher  than  case  III.  While  both  image  segments  have 
the  same  range  of  pixel  values,  the  "surface"  represented  in  case  I  fills 
much  more  space  than  that  of  case  III,  which  is  a  relatively  steep,  albeit 
smooth,  surface.  The  ambiguity  arises  when  case  11  is  compared  with  case 
III.  While  the  "shape"  of  case  11  is  much  rougher,  the  "range"  of  case  HI  is 
greater;  since  both  the  shape  and  the  range  of  an  image  segment  contribute 
to  its  roughness,  it  is  impossible  to  find  an  operator  that  assigns  a  fractal 
dimension  consistent  with  both  criteria. 


5.  Scaling  Pixel  Values 


We  may  try  to  avoid  this  ambiguity  by  assigning  a  fractal  dimension  based 
entirely  on  the  shape  of  the  image.  To  do  this,  we  must  scale  all  the  pixels 
in  the  window  to  fall  within  a  certain  range,  for  example  the  interval  [-10, 
10],  with  the  central  pixel  assigned  a  value  of  zero.  'line  results  of  scaling 
the  three  cases  in  figure  1  are  shown  in  figure  2. 

Once  the  image  segments  have  been  scaled,  the  computed  fractal  dimen¬ 
sion  will  be  based  entirely  on  the  shape  of  the  image,  not  on  the  range  of 
pixel  intensities,  so  cases  I  and  II  will  have  identical  fractal  dimensions. 
Similarly,  as  both  scaled  cases  I  and  II  are  rougher  than  case  III,  they 
should  have  a  higher  fractal  dimension  than  the  smoother  surface. 

Since  it  may  not  be  desirable  for  cases  I  and  II  to  have  the  same  fractal  di¬ 
mension  (based  on  the  greater  pixel  range  of  case  I)  once  the  fractal  dimen¬ 
sions  are  computed  for  the  scaled  images,  they  could  be  rescaled  by  a 
weighting  factor  based  on  the  range  of  the  nonscaled  window.  This  would 
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Figure  2.  Image  segments  in  figure  1  scaled  to  interval  [-10, 10]. 
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cause  case  I  to  again  have  a  higher  dimension  than  case  II  (as  we  would 
expect),  but  it  may  allow  a  smooth  surface  with  a  greater  range  to  have  a 
higher  fractal  dimension  than  a  rough  surface  over  a  smaller  range. 

The  decision  to  scale  the  pixel  values  or  to  rescale  the  computed  fractal  di¬ 
mensions  must  be  based  on  the  image  source  itself,  and  some  balance  be¬ 
tween  these  may  be  appropriate.  For  our  application,  where  all  pictures 
had  approximately  the  same  d5mamic  range,  we  found  it  was  useful  and 
consistent  to  scale  every  local  window  to  the  same  range  and  to  base  the 
computed  local  fractal  dimension  entirely  on  the  shape  of  the  surface. 


6.  Linear  Regression  Procedure 

In  the  mass  scaling  method  of  computing  local  fractal  dimension,  cubes  of 
increasing  size  are  considered  about  the  central  pixel,  and  for  each  size 
cube,  the  number  of  (scaled)  pixels  it  contains  are  coimted  (Moghaddam, 
1991,  and  Voss,  1986).  For  example,  if  we  were  to  consider  a  9x9  window 
about  the  central  pixel,  we  would  first  scale  all  the  pixel  values  so  that  the 
central  pixel  has  a  value  of  zero,  and  all  others  fall  within  the  interval  [-4.5, 
4.5].  We  could  then  consider  "cubes"  about  the  central  pixel  of  sizes  1, 3, 5, 
7,  and  9,  and  we  could  count  the  number  of  scaled  pixels  contained  in  each 
cube.  The  cube  of  side  1  would  necessarily  contain  the  central  pixel  (since 
it  was  scaled  to  zero).  The  cube  of  side  3  would  contain  the  central  pixel 
and  all  pixels  in  the  3x3  grid  about  the  central  pixel  with  values  in  the  in¬ 
terval  [-1.5, 1.5].  Thus,  for  n  =  1, 3, 5, 7,  and  9,  we  could  count  N(n),  which 
is  the'  number  of  pixels  in  the  nxn  grid  about  the  central  pixel  that  fall  in 
the  interval  [-n/2,  n/2].  Based  on  the  scaling,  N(l)  =  1  and  N(9)  =  81,  and 
N{n)  is  a  nondecreasing  function. 

The  temptation  at  this  point  is  to  borrow  from  the  technique  used  to  com¬ 
pute  global  fractal  dimension  (Falconer,  1991)  and  take  the  natural  logs  of 
the  n  and  N(n)  values,  plot  them  against  each  other,  use  linear  regression 
to  find  the  best-fit  line,  and  take  the  slope  of  that  line  as  the  local  fractal  di¬ 
mension  of  the  central  pixel.  Again,  we  use  an  example  to  show  why  the 
mass  scaling  method  for  computing  local  fractal  dimension  so  often  fails. 

Consider  two  surfaces,  one  smooth  and  one  as  rough  as  possible.  For  the 
smooth  surface,  N(l)  =  1,  N(3)  =  9,  N(5)  =  25,  N(7)  =  49,  and  N(9)  =  81,  in 
that  every  possible  pixel  is  contained  within  every  cube  of  increasing  size. 
The  most  pathologically  rough  surface  is  one  where  no  pixels  (except  the 
central  one)  are  contained  within  any  of  the  cubes  of  size  less  than  nine,  so 
that  N(l)  =  N(3)  =  N(5)  =  N(7)  =  1,  and  N(9)  =  81.  Table  1  shows  a  chart  of 
the  n  and  N{n)  values  for  these  two  surfaces,  and  figure  3  shows  the  regres¬ 
sion  lines  for  the  natural  logs  of  n  and  N{n)  for  the  smooth  and  rough  sur¬ 
faces,  respectively. 

For  a  smooth  surface,  the  equation  of  the  best-fit  line  for  ln(n)  and  ln[N(n)] 
is  y  =  2x,  and  the  assigned  local  fractal  dimension  is  2.  However,  for  the 
rough  surface,  where  we  would  expect  a  fractal  dimension  near  3,  we  find 
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Table  1.  Box  counting 
for  scaled  pixel 
values. 


Figure  3.  Plots  of  the 
best-fit  lines  for  the 
smooth  and  rough 
surfaces.  The 
abscissae  are  ln(n) 
and  the  ordinates  are 
ln[N(«)]. 
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the  best-fit  line  is  1/  =  1.765a:  -  0.831,  and  the  assigned  fractal  dimension  is 
1.765.  Clearly,  counting  boxes,  taking  logs,  and  performing  linear  regres¬ 
sion  does  not  suffice  as  a  local  fractal  dimension  operator,  regardless  of  the 
speed  of  implementation. 

To  solve  this  problem  with  linear  regression,  it  is  possible  to  force  the  re¬ 
gression  line  through  the  point  (ln(9),  ln(81)),  but  while  this  will  assign  a 
higher  fractal  dimension  for  rougher  surfaces,  it  sets  a  theoretical  upper 
bound  on  the  computed  fractal  dimension  of  2.82  for  a  9x9  window.  (In¬ 
creasing  the  window  size  increases  this  upper  bound,  possibly  beyond 
three,  and  therefore  this  is  an  artificial  "fix"  that  we  found  to  be  insuffi¬ 
cient  for  our  application.)  Instead,  we  borrow  a  procedure  from  the  box- 
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counting  method  used  to  compute  global  fractal  dimension,  and  propose 
this  as  a  new  method  of  computing  local  fractal  dimension. 


7.  A  Local  Fractal  Dimension  Operator 

In  order  to  compute  the  local  fractal  dimension  of  a  single  pixel,  a  window 
of  surroimding  pixels  must  be  selected.  While  the  literature  differs  on  the 
optimal  window  size,  ranging  from  9x9  (Moghaddam,  1991)  to  41x41 
(Keller,  1989),  we  must  achieve  a  balance  among  getting  enough  informa¬ 
tion,  having  the  procedure  be  computationally  tenable,  and  getting  so 
much  information  from  a  large  window  that  the  local  characteristic  of  the 
operator  is  lost.  This  optimal  balance  is  dependent  on  the  image  source, 
scale,  and  desired  speed  of  processing,  and  must  be  determined  for  each 
application.  For  the  algorithm  we  employ  for  our  application,  a  16x16  win¬ 
dow  about  four  central  pixels  was  found  to  be  appropriate. 

Since  we  use  a  window  of  even  size,  we  actually  compute  the  local  fractal 
dimension  for  a  spot  between  four  adjacent  pixels;  but  since  our  target  ob¬ 
jects  will  have  a  size  greater  than  one  pixel,  this  does  not  create  a  problem 
with  their  identification.  We  begin  by  scaling  to  zero  the  average  of  our 
four  central  pixels,  and  linearly  scaling  all  surrounding  pixel  values  in  the 
16x16  window  to  the  interval  [-8, 8].  We  then  compute  the  local  fractal  di¬ 
mension  of  the  center  point  using  a  method  similar  to  the  "box-counting" 
technique  for  computing  global  fractal  dimension,  and  different  from  the 
mass  scaling  method  discussed  in  section  6. 

The  center  of  the  16x16  window,  which  has  value  zero,  is  the  center  of  our 
domain  (i.e.,  the  plane  containing  the  pixel  positions).  The  range  of  scaled 
pixel  values  is  [-8,  8],  which  corresponds  to  heights  above  or  below  the 
domain.  This  cube  constitutes  the  space  used  for  computing  the  local 
fractal  dimension.  We  consider  cubes  of  size  1, 1/2, 1/4,  and  1/8,  and 
count  boxes  to  obtain  N(l),  N(1  /2),  etc.,  where  N(n)  is  the  number  of  boxes 
in  size  n  containing  one  or  more  scaled  pixels. 

Since  a  cube  of  size  1  covers  the  entire  space  (and  thus  all  scaled  pixel  val¬ 
ues),  N(l)  =  1.  For  this  method,  we  next  divide  the  space  into  eight  octants 
and  count  each  octant  that  contains  a  scaled  pixel  value  to  obtain  N(l/2). 
However,  to  proceed,  we  must  modify  the  box-counting  method  for  the 
reason  shown  m  the  following  examples. 

We  demonstrate  the  modified  box-counting  technique  with  a  two-dimen¬ 
sional  example  of  determining  N(l/8).  Consider  the  two  image  surface 
segments  in  figure  4,  which  have  been  scaled  so  that  all  pixel  values  fall  in 
the  interval  [-8,  8].  This  can  be  thought  of  as  a  slice  of  our  space  taken 
along  one  row  of  pixel  positions  in  the  domain.  We  have  overlaid  the  seg¬ 
ments  with  8x8  grids,  and  have  indicated  the  scaled  pixels  with  points. 

In  figure  5  we  cover  these  image  surfaces  with  boxes,  as  is  done  in  the  tra¬ 
ditional  box-coimting  method;  i.e.,  every  box  containing  a  pixel  is  counted 
to  compute  N(l/8).  Clearly,  while  the  image  in  case  I  is  more  "rough," 
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Figure  4.  Two 
segments  of  scaled 
image  windows. 
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Figure  5.  Box-covering 
of  scaled  image 
windows.  Ml/S)  =  9 
for  case  I,  while 
Ml/8)  =  16  for  case  II. 


Case  I 


Casell 


N(l/8)  is  greater  for  case  II,  which  leads  to  a  higher  computed  fractal  di¬ 
mension.  Although  the  surface  described  by  case  I  may  fill  more  of  the 
space,  the  discretization  of  the  surface  shows  pixel  values  as  points  with 
no  "fill"  between  them.  Additionally,  since  we  are  computing  the  fractal 
dimension  of  a  point  that  has  been  scaled  to  zero,  pixel  values  of  a  greater 
magnitude  should  indicate  a  higher  local  fractal  dimension. 


We  therefore  modify  our  box-coimting  technique  as  shown  in  figure  6, 
where  we  include  in  our  count  all  the  boxes  between  the  pixels  and  the 
domain.  This  technique  of  "stacking"  boxes  from  the  plane  of  the  domain 
not  only  measures  the  deviation  from  the  central  location,  but  it  also  mea¬ 
sures  the  amount  of  space  filled  by  the  surface  in  the  local  window. 

We  use  this  modified  box-covmting  technique  ("stacking"  boxes  from  the 
plane  of  the  domain)  to  compute  both  N(l/4)  and  N(l/8). 

Since  we  know  the  abscissae  of  the  regression  line  (ln(l),  ln(2),  ln(4),  and 
ln(8)),  we  express  the  fractal  dimension  slope  (FD)  of  this  line  entirely  m 
terms  of  N(n)  in  the  following  manner: 


FD  = 


31n  N 
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Figure  6.  Modified 
box-counting  by 
"stacking"  boxes. 
NowN(l/8)  =  21for 
case  I,  and 

Ml/8)  =  19  for  case  II. 


Case  I 


Additioirally,  N(l/2)  e  [4,  8],  N(l/4)  e  [17,  64],  and  N(l/8)  e  [67,  512]. 
(Theoretically,  the  lower  bounds  on  N(l/4)  and  N(l/8)  are  16  and  64,  re¬ 
spectively,  but  since  we  scaled  the  pixel  values  to  [-8,  8]  and  at  least  one 
pixel  takes  on  a  scaled  value  of  an  endpoint,  this  forces  these  lower  bounds 
to  be  slightly  higher  for  flat  surfaces.)  We  may  include  a  step  in  the  algo¬ 
rithm  to  identify  "completely"  flat  surfaces,  to  which  we  assign  a  local 
fractal  dimension  of  2.0.  This  method  clearly  places  our  computed  local 
fractal  dimension  in  the  interval  [2,  3].  (Naturally,  for  a  computer  imple¬ 
mentation  where  speed  is  important,  the  division  by  10  ln(2)  is  omitted  in 
the  interest  of  computing  time.) 


8.  Implementation  and  Modifications 

To  make  the  final  characterization  of  an  object,  we  start  with  a  bank  of  test 
images  to  which  we  apply  our  algorithm.  We  then  create  a  histogram  of 
the  local  fractal-dimension  values  within  the  objects  of  interest,  from 
which  we  create  a  look-up  table  based  on  an  appropriate  range.  We  search 
our  processed  images  for  these  values  to  detect  possible  objects  of  interest. 
Methods  exist  that  significantly  reduce  the  search  time  (Anderson,  1993), 
and  these  may  also  be  applied  to  reduce  the  number  of  pixels  for  which  we 
actually  compute  the  fractal  dimension;  however,  these  are  not  critical  to 
the  presentation  of  our  algorithm. 

While  our  algorithm  is  tailored  to  a  specific  application,  we  have  presented 
the  procedure  as  generally  as  possible  so  that  it  might  be  tailored  for  other 
applications.  This  is  best  done  through  experimentation,  the  proof  of 
which  is  the  final  result.  In  our  case,  we  want  the  fractal  dimension  of  a 
manmade  object  to  be  distinguishable  from  the  desert  sand  (and  perhaps 
even  sensitive  enough  to  distinguish  among  classes  of  vehicles);  but  the 
range  of  our  look-up  tables  must  allow  us  to  achieve  an  appropriate  bal¬ 
ance  between  false  positive  and  false  negative  readings.  Since  local  fractal 
dimension  is  not  unique,  false  positive  detections  are  always  possible,  but 
too  narrow  a  window  in  our  look-up  tables  may  lead  to  false  negative 
readings,  which  are  unacceptable  for  our  application. 
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A  16x16  window  size  is  appropriate  for  our  application,  based  on  the  im¬ 
age  source  and  a  specific  resolution,  but  this  can  certainly  be  modified  as 
the  source  changes.  For  example,  as  the  altitude  of  a  passive  image  source 
increases,  the  objects  of  interest  contain  fewer  pixels,  and  it  may  be  appro¬ 
priate  to  decrease  the  size  of  the  local  window.  Additionally,  since  this  lo¬ 
cal  fractal  dimension  is  not  scale-invariant,  we  need  to  recalibrate  our  look¬ 
up  tables.  Clearly,  for  other  applications,  the  optimal  window  size  can  be 
determined  through  experimentation. 

We  scale  each  window  to  the  interval  [-8, 8]  so  that  the  local  fractal  dimen¬ 
sion  is  based  entirely  on  the  "shape"  of  the  object,  and  we  do  not  re-scale 
these  computed  values  based  on  the  range  of  original  pixel  values.  For 
black-and-white  photography,  which  is  sensitive  to  different  light  condi¬ 
tions,  this  may  not  be  appropriate.  Instead,  the  computed  values  can  be 
scaled  by  an  appropriate  factor  based  on  the  dynamic  range  of  pixel  values 
as  mentioned  above. 

We  can  further  refine  our  procedure  to  make  it  even  more  sensitive  to 
slight  differences  in  pixel  values,  although  this  may  degrade  the  robust¬ 
ness  of  the  algorithm  for  some  purposes.  Since  strictly  counting  boxes  pro¬ 
duces  the  same  result  whether  the  greatest  pixel  value  in  our  square  is  6.1 
or  8.0,  we  may  increase  the  sensitivity  by  using  fractions  of  boxes.  For  ex¬ 
ample,  while  a  pixel  value  of  8  requires  four  boxes  of  side  2  to  cover  it,  a 
pixel  value  of  6.1  requires  only  3.05  boxes,  so  we  can  use  that  value  in  our 
algorithm.  To  ensure  that  our  computed  fractal  dimension  remains  be¬ 
tween  two  and  three,  however,  we  should  use  fractions  of  boxes  only 
when  more  than  one  box  is  required  to  cover  a  set  of  pixels.  This  modifica¬ 
tion  may  be  appropriate  for  applications  that  require  increased  sensitivity. 

If  we  are  concerned  with  actually  assigning  a  value  to  each  pixel  in  our 
image,  we  can  take  the  average  of  the  four  surrounding  computed  dimen¬ 
sions  as  the  value  we  assign,  but  for  most  applications  where  speed  is  im¬ 
portant,  this  step  is  unnecessary. 

A  disadvantage  of  this  procedure  is  that  it  leaves  a  band  eight  pixels  wide 
aroimd  the  border  of  the  image.  The  algorithm  can  be  modified  for  pixels 
near  the  border,  but  with  an  image  such  as  an  x-ray  where  the  area  of  con¬ 
cern  is  generally  in  the  center  of  the  image,  or  in  a  series  of  aerial  images 
with  overlap  between  adjacent  images,  we  need  not  worry  about  the 
border. 

This  procedure  produces  a  higher  fractal  dimension  near  the  edges  of  ob¬ 
jects,  and  while  our  goal  is  not  to  create  an  edge  detector,  this  feature  may 
be  useful  in  our  attempt  to  identify  objects.  As  the  fractal  dimension  of 
natural  features  varies  smoothly  across  an  image  (Super,  1990),  the  fact 
that  this  procedure  has  jumps  near  the  edges  of  objects  can  be  used  to  the 
advantage  of  the  analyst,  particularly  when  he  or  she  is  attempting  to  dis¬ 
tinguish  man-made  objects  from  a  natural  background. 
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